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Abstract. 

We present a trajectory-based semiclassical calculation of the full counting 
statistics of quantum transport through chaotic cavities, in the regime of many 
open channels. Our method to obtain the mth moment of the density of 
transmission eigenvalues requires two correlated sets of m classical trajectories, 
therefore generalizing previous works on conductance and shot noise. The 
semiclassical results agree, for all values of m, with the corresponding predictions 
from random matrix theory. 



PACS numbers: 05.45.Mt, 73.23.-b, 03.65.Nk 

1. Introduction 

Phase-coherent electron transport through ballistic quantum dots display a number of 
universal properties if the corresponding classical dynamics is chaotic [1]. These are 
well described by random matrix theory (RMT), in which the system's scattering 
matrix is assumed to be a random element of the appropriate ensemble [2], i.e. 
to be random unitary or unitary symmetric depending on whether time reversal 
symmetry is absent or present, respectively. RMT is therefore concerned with the 
average behaviour of a collection of different systems within a given universality class 
determined solely by the existing symmetry. On the other hand, properties of generic 
individual chaotic systems are expected to agree with the predictions of this theory 
provided a local energy average is considered, around a classically small but quantum 
mechanically large energy window. Rigorously deriving this connection between chaos 
and universality is one of the main challenges of the semiclassical trajectory-based 
approach to mesoscopic transport [3-5] . 

As in the case of spectral statistics of closed chaotic systems [6], the main 
ingredient from the classical dynamics is the existence of correlations between long 
trajectories. They organize themselves into families according to their action, and the 
elements of a family differ among themselves only by their behaviour in small regions 
(much smaller than their total length) in which some of them have crossings while 
others have anticrossings. The correlations induced by the existence of crossings is 
responsible for the emergence of universal quantum properties. This approach has 
been successful in reproducing RMT results for the average conductance [4,7], shot 
noise [8-10], time delay [11], and other properties [5]. All these calculations have a 
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natural perturbative structure in which 1/N, the inverse number of open quantum 
channels, plays the role of a small parameter. 

In this work we advance this line of investigations by obtaining, using the 
semiclassical approximation and classical correlations, the full counting statistics 
of chaotic cavities: the complete set of moments of the density of transmission 
eigenvalues. Physically, this encodes information about the statistics of the electric 
current through the system [12-14], viewed as a random time signal. The first two 
such moments are related to the average conductance and shot noise. Our result 
is restricted to leading order in 1/N, and to this extent we conclude that all linear 
statistical information contained in the RMT of quantum chaotic transport is also 
contained in the semiclassical approximation. 

In order to help put our results in perspective, let us draw an analogy with closed 
systems. By far the most popular quantity to be calculated in that case is the spectral 
form factor, the Fourier transform of the 2-point correlation function. The calculation 
requires taking into account contributions due to pairs of periodic orbits. By contrast, 
in the present work we consider two sets of trajectories with to elements each for 
all values of to. In a closed system this would correspond to obtaining all TO-point 
correlation functions. 



2. Counting statistics of chaotic transport 

Quantum transport is governed by the transmission matrix t, or equivalently the 
hermitian matrix tv. This matrix has a set {T l7 . . . ,T n } of n = mm{Ni, N 2 } non- 
zero transmission eigenvalues, where Ni and N 2 are the number of open channels 
in the incoming and outgoing leads, respectively. We consider a chaotic cavity 
with typical dwell time td, Lyapunov exponent A^ and linear size L. Together 
with the Fermi wavelength Xp, these last two quantities define the Ehrenfest time 
te = AT/ log(L/Air), roughly the time it takes for an initially localized wave packet 
to spread to the size of the system. For times much longer than te one can expect 
any initial state to become effectively equidistributed. 

The regime we are interested in is the semiclassical limit Ai? < I, when there 
are many open channels, Ni, N 2 1. However, we must approach this limit in such 
a way that t d , which is a classical time scale, satisfies t d ^> t E - Since te grows only 
very slowly as A^ — > 0, one may think of a suitable simultaneous shrinking of the 
width of the leads. This is the regime in which universality due to chaos is expected 
(for studies considering the situation when te > td, see [15] and references therein), 
and RMT predicts that the transmission eigenvalues are distributed in the interval 
[1 — 4£, 1] with average density given by [2, 16] 



P (T)=" ^-1. (1) 
Here N — N\ + N 2 is the total number of channels and the variable £ is defined as 

NN 2 

f=-^- ( 2 ) 

This result is valid to leading order in iV _1 , and is the same for all universality classes. 
The presence or absence of time-reversal symmetry is only felt in higher-order terms, 
sometimes called 'weak-localization' corrections. 
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The function p(T) can be characterized by its moments, 

M m = J p(T)T m dT = (Tr[(ttt)-]}, (3) 

where the brackets denote an average over the corresponding random matrix ensemble. 
The first two moments are related to the conductance (oc Mi) and to the shot noise 
(oc Mi — M2). For general m, the RMT moments were calculated explicitly [17], 

m-l , _\\ 

M m = N^[ p )(-l)%e P + 0(l), A^ao, (4) 



p=0 



where c p = ( 2 ^) /(p+ 1) are the Catalan numbers. Alternatively, they can be encoded 
in the generating function [18] 



2s \ V 1 - s 



This is the function we shall obtain scmiclassically. 



3. The semiclassical approximation 

A semiclassical approximation to the ti element of the transmission matrix is available 
from the corresponding theory for the Green's function [3] . It is given as a sum over 
all trajectories connecting incoming channel i to outgoing channel o, 

Uo = J2 (6) 

where as usual S 1 is the action of trajectory 7 and A 1 is related to its stability. The 
semiclassical expression for the moments contains 2m sums over classical trajectories, 

m 

^ c =EEn^^(^- s ' )/ft )- (7) 

i.o 7,<r j = \ 

Here i — . . . , i m } and o = {o\, . . . , o m } arc sets of m incoming and outgoing 
channels, respectively. 

We are denoting by 7 and <r two sets of m scattering trajectories: 7j goes from 
ij to Oj while o~j goes from ij+i to Oj (i m +i = h)- The quantities S-y = Y]j S 7j and 
S„- = J2j are the total actions of these sets. The average in (7) is taken around a 
classically small energy window, small enough to keep the classical dynamics and the 
amplitudes A 1 essentially unchanged. The phase factor on the other hand is rapidly 
oscillating as h — > 0, and this selects from the sum only correlated trajectories, with 
total action difference of order h. We restrict ourselves to the leading order in A^ 1 . 

In what follows we identify the classical trajectories possessing the required 
correlations and therefore contributing to (7). These are sets of trajectories involving 
encounters, examples of which are shown in Fig.l. A pair (7,0") will contribute to (7) 
only if there are places where 7-trajectories come very close to each other, forming 
an "encounter", and the corresponding cr-trajectories are obtained by a process of 
'reconnection' at the encounters, such that cfj initially runs closely to "fj+i and ends up 
running closely to jj . An encounter involving i trajectories is called an ^-encounter. In 
between encounters the trajectories follow arcs, along which the two sets are practically 
indistinguishable. The duration of a typical encounter is of the order of the Ehrenfest 
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Figure 1. Schematic examples of correlated sets of classical trajectories 
contributing to M| c . Each trajectory jj goes from incoming channel ij to 
outgoing channel oj, and is represented by a solid line. Trajectories Cj, which 
go from ij+i to Oj, arc in dashed lines. The circles mark the encounters, where 
trajectories switch partners (see text). 



O 6 O 
02 13 3 

Figure 2. The trees that correspond to Figure 1. The edges are the common 
arcs, the empty circles are the lead channels, and the shaded circles ("nodes") 
are the encounters. The edge emerging from i± is the root, and the channels are 
leaves. Note that the leaves are ordered ii , oi , ii , . . . , o m , starting with the root 
and going anti-clockwise. 



time te , much smaller than the typical duration of an arc which is proportional to the 
mean dwell time Trj. The action difference between the two sets of trajectories comes 
almost entirely from the vicinity of the encounters, and thus becomes small in the 
semiclassical limit. This theory has been discussed in several previous semiclassical 
calculations [4,5,7-11]. 

To perform the calculation we must construct all possible sets 7, and to this 
end we represent each set by a diagram containing its 'backbone' morphology of arcs 
and encounters. The complicated arcs of the actual trajectories are represented by 
straight edges; the encounters are represented by vertices of even degree (^-encounter 
has degree 2t) and the lead channels are represented by vertices of degree 1. The 
former vertices will be called nodes and denoted by shaded circles, see Fig. 3. The 
vertices of degree 1 will be called leaves and denoted by empty circles. We will see 
that our diagrams happen to be of a special kind, namely rooted planar trees. 

To each of these diagrams we associate a vector v = («2,f3, ■ • where is 
the number of ^-nodes (or ^-encounters). For example, one of the graphs in Fig. 2 
has characteristic v = (3), and the other has characteristic v = (1, 1). If a diagram 
has characteristic v, it contains V(v) = J^t v l nodes, while the number of edges is 
L(v) = m + J2i^ v i- Simple rules have been established to 'read off' the contribution 



(a) 



O ii 



(b) 



O ii 
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of a given pair (7, er) to M^. Their derivation involves setting local coordinates at the 
encounters using Poincare sections and invoking ergodicity to write down a probability 
density for encounters leading to an action difference AS. Then integrating over AS, 
and the possible duration of the arcs and applying a sum rule of Richter and Sieber [4] . 
We do not repeat this procedure, which has been reviewed in detail in [5]. The result is 
that each arc contributes l/N, while each encounter contributes —N. In the following 
sections we will use these rules together with an enumeration of the contributing 
diagrams to recover the random matrix prediction (4) semiclassically. 

4. No coinciding channels 

Let us initially assume that all channels are distinct, and do some simple power- 
counting with the channel number. There is a total of = N™N™ possibilities 
for distributing the incoming and outgoing channels. Suppose we have found a pair 
(7, er) represented by a graph with characteristic v. The edges produce a factor of 
\/N L{ - v \ while the nodes produce (~N) v(v \ We thus have N^NJ 1 /N L ^- V{ - V \ We 
want our result to be of leading order, so we must maximize V(v) — L(v), where 

V = V + 2m is the total number of vertices (including encounters and channels) . The 
quantity V(v) — L(v) is the negative of the Euler characteristic of the diagram, and 
it is well known that its maximal value is 1. Our moments therefore scale linearly 
with the number of channels, in agreement with (4). Moreover, it is also known that 

V — L = 1 if and only if the diagram is a tree. 

As the root of the tree we choose the edge containing i\. The defining feature 
of the diagram is the existence of a traversal i\ — > o\ — > ii — > . . . — > o m — > i\, such 
that each edge is traversed exactly twice: once in each direction. This implies that all 
leaves below any given node ( "below" in the sense of the natural ordering with the root 
being on the top) are consecutive with respect to the above traversal. Inductively, one 
can conclude that the branches of the diagram can be arranged in such a way that the 
leaves are ordered ii, Oi,Z2, . . . , o m , starting with the root and going anti-clockwise. 
Conversely, any tree with Vj nodes of degree 2j and the leaves marked ii, 0\, 12, ■ ■ ■ , o m 
anti-clockwise represents a diagram. A trajectory 71, for example, can be read off a 
tree by going from i\ to o\ along the shortest path. Thus we have established that the 
diagrams contributing to the leading order are in a one-to-one correspondence with 
■planar rooted trees. The term "planar" refers to the fact that the tree is defined up to 
an orientation-preserving homeomorphism of the plane onto itself, i.e. swapping the 
branches (generally) changes the tree. 

The total contribution of such a tree to is simply (— l) y ( v ) N£ m . A tree with 
characteristic v has L(v) — V(v) + 1 leaves so we denote by d(v) — (L(v) — V(v) + l)/2 
the number of incoming/outgoing channels associated with trees of characteristic v. If 
A/"(v) denotes the number of trees characterized by v, then their combined contribution 
is N£ m C m with 

C m = J2 AA(v)(-l)^>. (8) 

v: d(y)—m 

The function N(y) was studied by Tutte [19], who found an explicit formula for it, 
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However, for the sake of introducing the machinery needed in the following section, 
we take a different route and compute the above sum by making use of generating 
functions. The main idea is to notice that if 

f(x 2 ,X 3 ,...) )^A'(v:,V (10) 

V 

is the generating function of J\f(v) then #(r), obtained by setting x n+ i = — r n , is the 
generating function of the numbers C m , 

oo 

g(r) = f(-r, -r\ . . .) - ^ AA(v)(-l) v ( v V d W = ]T C m r m ~\ (11) 

v m— 1 

A tree characterized by v contains subtrees emerging from the top (n + l)-node, 
which may be characterized by their own node vectors Vi, . . . , v 2 „+i, numbering left 
to right. One can therefore count the number of trees characterized by v by counting 
all possible subtrees. This implies a recursion relation for A/"(v) (see Appendix A for 
details), 

^( V )^E E 11^)5^-^, (12) 

n>lvi-v 2 „+i j=l 

where w = J^j w j an d e « has 1 in n-th entry and zero everywhere else. Substituting 
the recursion relation in Eq. (10) we see that / satisfies / = 1 + X2.f 3 + x%f 5 + . . .. 
Correspondingly g = 1 — rg 3 — r 2 g 5 — . . .. Summing the geometric series we arrive at 
,g = l — rg 2 , and thus 

-1 + VTT4F 

ffW = • ( } 

When compared with the generating function for the Catalan numbers c m this gives 

C m = (-l)™- 1 ^-!, (14) 
in agreement with Eq. (4). 



5. Coinciding channels 

The previous calculations assumed all channels to be different. If it happens that 
ij = ij+\, then jj and 7j+i enter the cavity from the same channel. It is useful to 
view the corresponding diagram as arising from a more general one (with ij ^ ij+i), 
in the limit when a particular encounter happens closer and closer to one of the leads. 
After taking this formal limit, all trajectories previously involved in the encounter now 
enter from the same channel, and the encounter has disappeared. This is illustrated 
in Fig. 3 for a simple example with m — 3. We are actually neglecting any spatial 
dimension the encounter may have and treating it as point-like. This is justified 
insofar as we consider te <C td , since the length and width of encounters are typically 
proportional to the Ehrenfest time. 

In the situation just described we say that the node "touches" the lead, which may 
be incoming or outgoing. In any case, it no longer gives the usual —N contribution. 
We now have to count in how many ways a given one of our trees can have nodes 
touching leads. Consider for example the tree in Fig. 2(a). The lower node has two 
edges leading to outgoing leaves. This node could touch the outgoing lead, and hence 
these edges would actually vanish. The top node could touch the incoming lead, while 
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Figure 3. The diagram on the right, which has coinciding channels, is seen 
as a limiting case of the more general diagram on the left. The encounters are 
continuously brought up to the leads, and eventually disappear as all involved 
trajectories enter (or leave) from the same channel. 




Figure 4. Left: Schematic of trajectories contributing to M| c when iz = 14- The 
crossing 'touches' the incoming lead. Right: The corresponding tree. Compared 
to the tree in Fig. 2(b) the number of segments is reduced by two, and the number 
of crossings is reduced by one. 



the middle node cannot be made to touch any of the leads. In Fig. 4 we show the 
schematics and the corresponding tree of a contribution similar to Fig. 1(b) but with 
h = U- 

When an (n + l)-node touches an incoming lead, the number of edges is reduced 
by n+ 1. The number of nodes is reduced by 1, so the contribution gets multiplied by 
— N n . Because the number of channels involved in the sum is also reduced, there is a 
factor of N^ n . If we define z\ — N/N\, making an (n + l)-node touch the incoming 
lead produces a factor — z™ multiplying the contribution of the tree. For a given tree 
T, denote by qi iTl (T) the number of (n + l)-nodes that can be made to touch the 
incoming lead. There are ways to have k out of the qi tTl such nodes actually 

touch the lead. Therefore, the basic contribution must be multiplied by 

E(-^) fc ( 9l : n )=(i-^) ai -- (is) 

fe=0 ^ ' 

It is easy to see that each and every diagram with coinciding channels can be derived 
from one and only one "parent" diagram, which does not have any coinciding channels. 
The conclusion is that we may consider only these "parent" diagrams, provided we 
multiply each one of them by (1 — z") 9l '"(l — z^) 92 -", where the subscript 2 refers 
to the outgoing lead. This accounts for all possible contributions in which channels 
coincide. 



Full counting statistics of chaotic cavities from classical action correlations 



8 



The above arguments amount to saying that, similarly to (8), the complete 
semiclassical moments are given by 

M-(0 = JVr ]T (-l) y ^P(z 1 ,z 2 ;v), (16) 

v: d(v)—m 

where the sum is over all characteristic vectors v. Instead of including only the number 
7V(v), as in (8), we must instead use 

p( Zl , Z2 ;v)= n( 1 -^r" ,(T) (i-^r- (T) ; (17) 

TeT v n>l 

where the sum is over T v the set of all trees characterized by v. Notice that 
z\z 2 = z\ + z 2 = Now we define the generating function of P, 

F(z 1 ,z 2 ;x 2 ,X3, . . .) = ^2P( Zl ,z 2 ;v)x v 2 2 x v 3 3 (18) 

V 

and, by defining G sc (r) = N£F(zi, z 2 ; —r, — r 2 , . . .), we obtain the analogue of (11) as 

G sc (r)=N^^ Lrm ~ 1 =J2 MS ^h) ■ (19) 

Unfortunately, it is not straightforward to write a recursion for P. This is because, 
unlike any other node in the tree, the top node can be made to touch both leads 
(although not at the same time). To circumvent this problem, we define an auxiliary 
function, 

P( Zl ,z 2 ;v)= n( 1 -^) 9 ^ (T) (l-^) g2 - (T) ; (20) 

TeT v n>l 

where q[ „ is defined as q ijn but excluding the top node. For P there is a natural 
recursion relation, analogous to (12), explained in the appendix. It is given by 

?{zi,z 2 ;\) = ^1^2; vi) 

n>l Vi---V2n + 1 

(21) 

where again we use w = X^o 1 v i ano - u = Sj=o v 2j+i- The sum is over the valency 
of the top node and the characteristic vectors of the 2n + 1 subtrees that emanate 
from it. The factor z 2 <5 Uj o includes the contribution due to the top node touching 
the outgoing lead which is only possible when all odd Vj's vanish. The function P is 
useful because it is related to P according to 

n 

P(z u z 2 ;v) = Y J P(zi,z 2 ; Vl )l[P( Z2,zi;v 2 j) 

n>l Vi-V2 n +i j = l (22) 

xP(z 1 ,z 2 ; v 2 j+i)(l - z'i 5 w ,u - z 2 <5 u ,o) <5 W;V _ e „ +1 • 

Which now includes the contribution due to the top node touching the incoming lead 
when all the even Vj 's vanish. The details are also left to the appendix. 

We now denote by f(z\,z 2 ; x) the generating function of P. The recursion relation 
(21) implies 

f(z u z 2 ; x) = 1 + x n+1 (f n+1 (z u z 2 ; x) - z%)f n (z 2 , Zl ; x). (23) 

n>l 
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h= ^^' ' VVS -A^S'-TS U (26) 



We again identify x n+ i = — r n and write g(zi,z 2 ;r) — f{z\, z 2 ; —r, — r 2 , . . .). Making 
use of the geometric series, we can reduce (23) to 

g(zi, z 2 ; r) = 1 + r(z 2 - l)h(z 1 , z 2 ;r), (24) 

where we have defined a function which is symmetric with respect to the variables 
zi, z 2 , namely h(zi, z 2 ; r) — g(zi, z 2 ; r)g(z 2 , z\\ r). Substituting (24) back into h leads 
to the algebraic equation 

£h = £(1 - rhf + rh, (25) 

which can be solved to give 

£ + 2£r - r - ^/g - r)(4£r + £ - r) 
2£r 2 

Having this solution we can compute g(z\,z 2 ;r). Because of (22), the function 
/(zi,z 2 ;x) is related to F(z 1} z 2 ;x) simply by F = f - J2 n >i x n+iZif n+1 - This 
implies that 

G-(r) = m^lll (27 ) 
1 - rzig{zi,z 2 ;r) 

Let us multiply and divide the above expression by (1 — rz 2 g(z 2 , z\\ r)). We may 
then use that 

f - r 

(l-rz 1 g(z 1 ,z 2 ;r)){l-rz 2 g(z 2 ,z 1 ;r)) = — — , (28) 

by virtue of (24) and the definition of z\p. On the other hand, using (24) again 

g{zi,z 2 ; r)(l - rz 2 g(z 2 , z\;r)) = \- rh. (29) 
Finally, taking the explicit formula (26) into account and writing s = r/£ leads to 



E 




Mm(£)s m = - | x j 1 + ^ - 1 ] , (30) 

rn— 1 

which is identical to (5). Therefore, all semiclassical moments are indeed equal to 
the corresponding random matrix theory predictions. From the semiclassical point of 
view, the fact that the leading order result is the same with or without time-reversal 
symmetry stems from the fact that, being trees, the contributing diagrams contain no 
loops (in contrast to Richtcr-Sieber pairs [4], for example). 

Obtaining higher-order terms in iV _1 , which are necessary to describe experiments 
in which a relatively small number of channels arc involved, would require the 
incorporation of more general sets of trajectories, for which the corresponding 
diagrams are no longer trees. Computing the full perturbative series seems to be a 
very demanding task, especially for systems with time-reversal symmetry. We should 
mention that some exact RMT results (valid for arbitrary channel numbers) have 
appeared recently [20,21]. In particular, for broken time reversal symmetry the average 
values of TrT m , (TrT 1 ™) 2 and (TrT) m , where T = tt\ were obtained for general m 
by one of the present authors [22]. 
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Appendix A. Tree recursion relations 

Let T be a tree and T v the set of trees with characteristic v. We first establish the 
recurrence relation (12) for 7V(v) = \%r\, the number of trees characterized by the 
node vector v. To derive the relation we break the tree at the top (n + l)-node 
adjacent to the root. The top node has degree 2(n + 1) and splitting the tree T at 
this point the node becomes the root of 2n + 1 subtrees T\, . . . , T 2n+ i, characterized 
by vectors vi, . . . , v 2n+1 . Clearly v = X^=t v i + e « where e„ has 1 in its nth entry 
and zero elsewhere, representing the top node that was removed. Figure Al shows a 
tree with characteristic (4). This tree splits the at the top node, degree four (n = 1), 
into three subtrees each characterized by (1). In general the number of trees with top 
node degree 2(n+ 1) is given by the number of combinations of subtrees, rij=t A/"(vj) 

where w = X^=i v j = v — e «- Summing over all allowed valencies of the top node 
establishes the recursion relation, 

^(v) = £ E LI ^(v,-)<W-e„ +1 • (A.l) 

n>l Vl--V2n + 1 j = l 




O] i 2 o 2 i 3 o 3 ; 4 o A i 5 o 5 

Figure Al. A tree with characteristic v = (4) separates at the top node into 
three subtrees characterized by vi = V2 = V3 = (1). 

The recursion relation (21) for P can be generated in a similar manner. To recap, 
P(z 1 ,z 2 ;v)= ^2 ]J(l-z™)«U( T )(i- z ™)^, m (T) ; (A2) 

TeT v m>l 

where q2, m (T) is the number of ways an (m+ l)-node can touch the outgoing lead and 
q[ m (T) the number of (m + l)-nodcs - excluding the top node adjacent to the root - 
that can touch the incoming lead. For our trees the root corresponds to the incoming 
lead and our definition of P excludes the possibility of touching the root. 

To establish the recursion relation we again consider breaking the tree into 
subtrees Ti, . . . , T 2n +i at the top (n + l)-node, numbering left to right. The number 
of ways an (m + l)-node of the tree can touch the incoming lead is 

n 

<l'i, m (T) = gi im (Ti) + ]T (<Z2, m (T 2j ) + <7i, m (T 2j+1 )) . (A.3) 
j=i 

This includes a change in the ordering of the leaves on the subtrees with even index. 
On subtrees with odd index the first (left most) leaf must be an outgoing lead, labeled 

0, while the first lead of an even numbered subtree is labeled with the incoming lead 

1, see for example figure Al. 
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To find the number of ways an (m+ l)-node can touch the outgoing lead, similarly 
one sums the q 2y m of the odd subtrees and adds the q[ m of the even subtrees where 
the o and i leaf labels must be exchanged. In addition if m = n the top node can 
contribute. The top node may be made to touch the outgoing lead if all the odd 
subtrees have characteristic 0, i.e. every odd edge ends in a leaf. Figure A2(a) shows 
a tree where the top node can touch the outgoing lead. Therefore, 

n 

(A.4) 

where u = X)j=o v 2j+i- Consequently P(z!,z 2 ;v) is expressed in terms of functions 




Figure A2. (a) A tree with characteristic v = (2) where the top node can be 
made to touch the outgoing lead, (b) A tree with characteristic v = (3) where 
the top node can be made to touch the incoming lead. 

P(-, •; Vj) generated by the subtrees, 
P(z 1 ,z 2 ;v) = -P( z i> z 2;vi) 

n>l Vi ■■■V 2 n + l 

(A.5) 

X Y[P( z 2,Zi;V2j)P(zi,Z 2 ;V2j+l){i ~ Z 2 ^u,o) #w,v-e„+i 

From P we recover P by calculating the contribution generated when the top 
node touches the incoming lead independently Recall the definition of P, 

P( Zl ,zr,v)= ]T n(l-^) 9l "" (T) (l-^ m r"" (T) ■ (A.6) 

TeT v m>l 

Comparing this with the definition of P (A. 2), we see that terms in the sum are 
identical except when the top (n+ l)-node of T can touch the incoming lead, in which 
case qi, n {T) = q'i_ n {T) + 1- If we let 1Z V C T v denote the set of trees that can touch 
the root the definition of P can be rewritten using P. 

P(z 1 ,z 2 ;v) = P(z 1 ,z 2 ;v)- £ z[ l [] (1 - zT) q '^ T) (1 - 4T 2 < ™ (T) , (A.7) 

Te1Z v m>l 

where n is determined by the degree of the node adjacent to the root. A tree is in 
1Z V if, counting left to right, the even branches of the top node all end in leaves, 
v 2 j = for j e {1, . . . , n} or equivalently w = u. Figure A2(b) shows a tree in 1Z V . In 
equation (A.7) the product inside the sum over 7Z V is the same as that appearing in 
the definition of P and we can again break the trees in 1Z V at the top node to write the 
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contribution in terms of subtrees. Combining this with the recursion relation (A. 5) 
for P, we obtain the following formula for P: 

n 

P(z 1 ,z 2 ;v) = P(zi,z 2 ;v 1 ) Y[P(z2,z 1 ;v 2j ) 

n>lvi-v 2 „+i j = l (A-8) 

XP(Z1,Z 2 ; V 2 j+l)(l - Z™(5 W , U - z% 5 u ,o) 6 w ,v-e n+1 
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